Effective Coulomb interactions in solids under pressure 
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Correlated materials are extremely sensitive to external stimuli, such as temperature or pressure. 
Describing the electronic properties of such systems often requires applying many-body techniques 
to effective low energy problems in the spirit of the Hubbard model, or extensions thereof. While the 
effect of pressure on structures and bands has been investigated extensively within density-functional 
based methods, the pressure dependence of electron-electron interactions has so far received little 
attention. As a step toward ab initio pressure studies for realistic systems within a setup of max- 
imally localized Wannier functions and the constrained random phase approximation, we examine 
in this paper the paradigmatic pressure dependence of Coulomb interactions. While compression 
commonly causes the "extension" of Wannier functions, and thus transfer elements, to grow, we 
find the - seemingly counter-intuitive - tendency that the bare Coulomb interaction increases under 
compression as well. We reconcile these behaviors by appealing to a semi-analytical tight-binding 
model. We moreover argue that, for this model, the requirement of maximal Wannier localization is 
equivalent to maximizing the Coulomb interaction matrix elements. We then apply the above first 
principles techniques to fee hydrogen under pressure. While we find our comprehension of the bare 
Coulomb interaction confirmed, the induced changes in screening strengths lead to an effective one 
band model with a Hubbard interaction that is non-monotonous under pressure. 

PACS numbers: 



I. INTRODUCTION 

The panoply of structural, orbital and spin degrees of 
freedom and the joint presence of important electronic 
Coulomb interactions cause correlated materials^ to be 
the realm neither of band-theory nor of model many-body 
physics each on their own. Therefore, in recent years, 
much ingenuity was invested into finding ways to merge 
the "realism" of the former with the accurate description 
of correlation effects of the latter. 

With the exception of the GW approximation^ 3 -^ to 
Hcdin's equations 2 ., the electronic structure approach 
GW+DMFT^ that combines GW with the dynamical 
mean-field theory (DMFT)-, and a recent proposal for 
a self-energy downfoldingi, this combining commonly 
amounts to extracting a low-energy one-particle Hamilto- 
nian from density-functional-theory-based methods, such 
as the local-density approximation (LDA)^^ supple- 
menting it with interaction terms [e.g. of the Hubbard- 
Hund (U and J) type] and to solving the resulting "real- 
istic model" with a chosen many-body technique. 

Hence, in this approach there are two intertwined 
preparative tasks that generate the many-body problem. 
The deducing of the low-energy one-particle part can be 
achieved e.g. by a tight binding fit of relevant bands, the 
downfoldin g i procedure within e.g. muffin-tin-based 
methods ) 11 ' 12 ! 13 or by the generation of (maximally lo- 
calize d 14 ! 15 ) Wannier functions^ 

The interaction parameters of realistic many-body 
models, in turn, are often chosen rather empirically than 
from a solid first principles basis, a fact that has caused 
many objections in the past. In particular, when track- 
ing properties as a function of an external parameter - 



pressure, in our case - the evolution of the interaction 
has mostly been discarded. 

This and the quest for going beyond mere qualitative 
results towards, eventually, the quantitative design of 
materials, highlights the need for accurate ways to de- 
termine all ingredients of realistic models in an ab initio 
fashion. Nowadays, the most popular methods for the 
computation of interaction matrix elements are the con- 
strained LDA technique^, and the constrained random 
phase approximation (cRPA)i£. 

A recent, and promising approach is the use of Wannier 
functions within the cRPA setup-^, which allows for a de- 
ducing of the one-particle and two-particle parts of the 
Hamiltonian on the same footing. Moreover, working in 
a localized Wannier type of basis is often a requirement 
of many-body approaches such as the DMFT— . As to 
the interaction matrix elements, the cRPA technique al- 
lows for a precise elimination of the screening channels of 
the chosen orbital subspace that constitutes the effective 
modelii. 

While these techniques have already been applied for 
the setting up of many-body models of some complex 
material a 19 i 20 ' 21 , a basic understanding of the effects of 
pressure on the Coulomb interactions within a Wannier 
setup is lacking. This is the aim of the current work. 

In a first part, we investigate a semi-analytical tight 
binding model of a one-band solid in one dimension, 
track transfer matrix elements, the bare (i.e. unscreened) 
Coulomb interaction, and the spread of the maximally lo- 
calized Wannier orbitals as a function of lattice spacing. 
Being able to access the decomposition of the maximally 
localized Wannier functions onto the tight-binding basis 
will allow to understand the surprising finding that under 
pressure the Coulomb interaction matrix clement aug- 
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merits, while, at the same time, transfer elements describ- 
ing the delocalization of electrons grow as well. As a more 
realistic example, we, second, apply the fully ab initio ap- 
proach of the cRPA within maximally localized Wannier 
functions 1 ^ to fee hydrogen, which is found to exhibit the 
explained generic behavior of the bare Coulomb interac- 
tion matrix elements. However, the partially screened 
Coulomb interaction - the Hubbard U of an effective low 
energy model for the half-filled Is orbital - actually shows 
a non-monotonous trend - a consequence of two opposing 
effects in screening processes. 



II. METHOD 

The method of using the Wannier orbital construction 
in conjunction with the cRPA technique has been pre- 
sented in detail in RefJ^. While not fully reviewing the 
approach, we will discuss some issues relevant for the un- 
derstanding of our results and introduce some notation. 

For the one-particle band-structure, a density func- 
tional calculation is performed. For the realistic case 
of fee hydrogen, we will employ the LDA^ in the full- 
potential linear muffin-tin orbital (FP-LMTO)^ 2 . realiza- 
tion. For obtaining the random phase approximation 
(RPA) polarizations we employ the code of Refi^ with 
the maximally localized Wannier extension of ReffiS, 
and construct an effective system for the isolated Is 
orbital. That is, we introduce the sub-Hilbert space 
7i off = span {iV'kis)} j spanned by the Is Kohn-Sham 
wave function. Since the aim is to work within a localized 
basis, the extraction of the low-energy part is done by a 
construction of Wannier functions 1 ^ for Tt° !1 , as described 
in Refsi 14 i 15 ' 19 . The corresponding effective interactions 
are then computed within the constrained RPAi£ formal- 
ism. This amounts to screening the matrix elements of 
the bare Coulomb interaction w(r, r') = 1/ |r — r'|, which 
in the Wannier basis are given by 
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R,R' 



(1) 



with a partial RPA polarization 

Pr = P ~ Ps 



(2) 



where P and P s are the polarizations of the full and the 
sub-Hilbert spaces, respectively. The latter P s , when us- 
ing the Kohn-Sham orbitals, can be expressed as 
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i.e. by transitions restricted to the effective sub-system, 
in our case Ti* tt = span {IV'kL)}- Within this notation, 
the strengths of screening channels arc influenced by two 
effects : the matrix elements [the overlap integrals of 
wave functions that occur when calculating matrix ele- 
ments of P, in analogy to Eq. JlJ] , and the energy differ- 
ences of the Kohn-Sham excitations, etn, that appear in 
the denominators. The virtue used for the constraining is 
the fact that the screening contributions are additive 1 ^. 
Indeed, when using the above decomposition P = P s +P ri 
the fully screened interaction W (of the GW formal- 
ism2iiL£) can be given in terms of the partially screened 
interaction for the effective model of the one-band or- 
bital, W r = v/(l-P r v), by W = W r /(1 - P s W r )±£. 
The Hubbard U of the Is sub-system is given by the 
on-site Wannier function matrix element of W r . 

The major observation in this context is that the con- 
struction of Wannier functions, and thus also the deter- 
mination of interaction matrix elements are not unique 1 ^. 
Indeed a unitary transformation applied to the periodic 
part of the wave functions, while preserving the Bloch 
functions, changes the Wannier orbitals. This gauge free- 
dom can be used to choose the Wannier basis that is most 
suitable for the final purpose. The aim of the low-energy 
system is to correct for local (Hubbard-type) interaction 
effects. To this end, there exist at least two proposals 
on how to choose an optimal Wannier basis set : in the 
maximally localized Wannier approac h 14 i 15 the extension 
("spread") of the Wannier functions is minimized : de- 
noting the Wannier states by kets, XR Q ( r ) — ( r l a R); this 
spread can be chosen as^ 4 - 
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The minimization of this O is of course only one of the 
possible options to fix the Wannier functions, yet a very 
natural one, since it can e.g. be show n 14 ' 24 that, in one di- 
mension, the resulting Wannier functions arc cigenfunc- 
tions of the subset projected position operator, an intu- 
itive criterion for real space localization^. 

In the second scheme, the screened local Coulomb in- 
teraction matrix element - the Hubbard U, as given 
by the on-site component of W r from above - is max- 
imized 1 ^^ to determine the Wannier functions. For the 
case of SrV03, it has been shown 1 ^ that both approaches 
yield very similar results. Indeed, appealing to e.g. the 
equation of the bare Coulomb interaction matrix element 
in the Wannier basis, Eq. it seems plausible that a 
greater localization of Wannier functions (smaller Q) re- 
sults in an increased interaction matrix element V. For 
the simple model in one dimension that we discuss in 
the following, we in fact motivate the equivalence of the 
maximally localized Wannier functions and the basis in 
which the (bare) Coulomb interaction matrix element is 
maximal. 

However, we stress that the intuitive correspondence 
between a stronger localization of Wannier functions, in 
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the sense of the spread CI, and a larger interaction ma- 
trix element in this basis does not hold in general. Here, 
one has to distinguish between the changes in the Wan- 
nier localization that are induced by a modification of 
the Wannier gauge from those caused by modifications 
of the lattice, caused e.g. by pressure. The fact that for 
the discussed model, both methods to fix the Wannier 
gauge arc equivalent, states that for a given pressure the 
maximally localized Wannier functions yield the maxi- 
mal Coulomb interaction matrix element. Yet, as we shall 
see, a pressure induced increase in the converged minimal 
Wannier spread is actually quite naturally concomitant 
with a greater bare interaction. 

III. RESULTS AND DISCUSSION 

External pressure, or structural changes in general, 
provides an impetus to alter not only the one-particle 
band-structure of a material, but also the Wannier func- 
tions. Therefore, when investigating the pressure depen- 
dence of effective, i.e. screened, interaction matrix ele- 
ments, one has to distinguish between influences of the 
former, which enter via a modification of screening chan- 
nels, and of the latter that not only affects the polariza- 
tion, but, on a more fundamental level, modifies already 
the bare Coulomb interaction matrix elements. The fact 
that Wannier functions of a solid are not eigenfunctions 
of the system, may result in a nonstraightforward evolu- 
tion when parameters such as external pressure change. 

A. One-band tight-binding model in one dimension 
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FIG. 1: (Color online) The discrete distribution A n for differ- 
ent lattice constants, a/ao, as a function of the atomic "dis- 
tance", n, in real space. The inset shows the dependence of 
Ao, i.e. the weight of the atomic function at the origin, on the 
lattice constant. 



t2 

with the eigenvalue — „ 5 . While the tight- binding 

approach with this orbital can in principle be used to 
approximate the Bloch eigenfunction for a Hamiltonian 
with any potential, its use is obviously most justified for 
a Kronig-Penney type of mode l 26 i 27 with a Dirac-comb 
potential. 

The Bloch function is written as : 

Mz) = A k J2e %kR x(x-R) (7) 

R 



As a first model system, we investigate a tight-binding 
parametrization of a one-band solid in one dimension. 
In that case the maximally localized Wannier function 
is naturally given by the Fourier transform of the Bloch 
functions when inversion symmetry is verifie d 14 ! 24 . More- 
over, since no higher energy orbitals are present, no par- 
tial screening can occur [in Eq. : P r = 0] and the 
Hubbard U equals the on-site matrix element of the bare 
Coulomb interaction, U = Vr.r. 



1. Bloch function in tight binding 

As building blocks of the tight binding basis functions 
we opt for a hydrogen-like Is orbital in one dimension 

1 



-|x|/o0 



(5) 



with the Bohr radius do- This orbital is a solution to 
the time-independent Schrodinger equation with a single 
binding delta potential : It represents an eigenstate of 
the "atomic" Hamilton operator 
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Here, the factors Ak assure the orthonormality of the 
Bloch function tpk, and is determined to be : 
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Sm cos(kma) 



(8) 



where, with the lattice constant a, and the integer m, ma 
denotes the distance to the mth neighboring site. Fur- 
ther, s m denotes the overlap integral between the atomic 
function at the origin and its mth neighbor, and is given 
by: 



(9) 



/ dxx*(x)x(x - ma) 
= fj [ |m|a ^ e -M»/qQ 

2. Wannier function 



In one dimension, the maximally localized Wannier 
function, tpji(x), is given by the usual Fourier trans- 
form of the Bloch function if inversion symmetry is vcri- 
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FIG. 2: (Color online) Wannier functions for different lattice 
spacings a/oo- The inset shows the deviation of the prob- 
ability distribution |-i/>(:r)| 2 from the atomic limit (note the 
different scales in the left and right inner panel). 



pressure is applied, and the lattice constant shrinks, fi- 
nite overlaps of the (non-orthogonal) hydrogen orbitals 
entails contributions from neighboring sites to mix in, 
and the distribution A„ broadens (see Figure [1]). Since 
^bn+i < 021, and |-A 2 | n |+i| > ^(H+i), the normaliza- 
tion of the Wannier function, Eq. (fl4|) . causes the coeffi- 
cient of the atomic orbital at the origin to become larger 
than 1 : Aq > 1. This results in a greater probability 
density |^r(x)| 2 around the site origins, x — R = 0, when 
pressure is applied. 

The corresponding Wannier functions of the above 
cases are shown in Figure [5] As anticipated from the 
above discussion, more weight is accumulated at the ori- 
gin : a harbinger for a larger on-site Coulomb interac- 
tion. On the other hand, contrary to the atomic limit, 
the tails of the Wannier functions extend over several lat- 
tice constants, before the exponential decay sets in. This 
behavior at larger distance points to an increase in the 
Wannier spread and the growing of transfer integrals. 



fiedii^i. Therewith 

1 f*/ a rlh 



= ^£f /a ^e-^- R U kX ( X -R') 



= Jj E A R-R'X(x - R') 

R' 

n 

where we defined 



(10) 

(11) 



T/Q ^l e -iHR-R>) Ak (12) 



— 7r/a 



2tt 



or, with R' = na, and the reference R — 0, 
A n _ 1 r /a dk ~ lkna 



C C J_ T/a 2tt v /l + 2E™ =1 s m cos(fcma) 

(13) 

Demanding J dx | (a?) | 2 = 1, the normalization con- 
stant C becomes : 



(14) 



The quantity A n is real for all n for symmetry reasons. 
FigurcQ]shows its behavior for different lattice constants : 
in the limit of large atomic separation (a 3> ao), the over- 
laps s m are negligible, and the Wannier function %jjo(x), 
Eq. (fTTjl . will equal the atomic orbital x( x )i Eq. ((5|). Thus 
the distribution A n picks up a single mode of the ar- 
ray, A n = 6 n o for the representative site "0" . When 



3. Wannier spread 

Since by symmetry, (r) = (x) = (R = 0\X\R = 0) = 0, 
the spread of the Wannier function [fl as given in Eq. Q] 
reduces to : 

(X 2 ) = (R = 0\X 2 \R = 0) = Jdxx 2 \iPr= (x)\ 2 



2 ( a , 2 2 

- + an 



+ 2 EE KA n+m e- ma 

n m>0 

/m \ 3 2 (ma) 2 

(y a ) 3^ + 4^ 



/ao 



■ao 



(15) 




FIG. 3: (Color online) Spread of the model as a function of 
lattice constant a/ao. The dotted line indicates the atomic 
limit : {X 2 } = a 2 /2. 
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While the first term is always positive, the second will 
be negative for odd m and positive for even neighbors 
(see the form of A n ). Given the decay behavior of A n , 
the second term will thus be negative in total, yet small. 
Indeed the major contribution to the spread comes from 
the first term, making it plausible that the spread, as 
defined by (X 2 ), increases with a more widely distributed 
A n , as shown in Figure [3J 

Although helpful for the understanding of the current 
model study, the spread is not too good a quantity for 
gaining quantitative insights from pressure studies of re- 
alistic systems, as we will discuss later. 



induced derealization increases the transfer integral, as 
expected. Figure [4] displays the hopping as a function of 
the lattice constant. 



5. Electron- electron interaction 

All of the above were concerned with the one-particle 
picture, i.e. the non- interacting Hamiltonian Hq. The 
on-site electron-electron interaction matrix clement Vo,o, 
that governs the two-particle term in the final many-body 
model, reads in the Wannier basis : 



4- Transfer integral 

For the transfer integral t, we need to explicitly specify 
the non-interacting Hamiltonian, Hq, and, for simplicity, 
we shall choose a Kronig-Penny-type mode l 26 i 27 with an 
ionic Dirac-comb potential, i.e. 



i p2 + 5>™ 

1 _ p 2_yJ^_ 5 r X _ la \ n 6) 



2m 



Then the nearest neighbor hopping t = (ipR=o\Ho\ipR=a) 
can be expressed as : 



2m e al t = ^ AnA 



m (2 S m _j_ i _ n 



(17) 



^ ' A n A m S ' 



\l— n\a/ao — \l— m— l\a/ao 



l^m+1 



As can be inferred from the dependence of the coefficients 
A n , Eq. (fT3|) . and the overlap s m , Eq. j9]), the pressure 



V 



Vo.o 



m e ao 



dx / dx' \ip (x)\ 2 v(x,x')\ip (x')\ (18) 



dx \ipo(x)\ l 



dx\"4>o{x)\ I dx 



,ivo(z+or 



(h.c.) 



(C.) 



when choosing a hardcore (h.c.) or Coulombic (C.) in- 
teraction : 



v(x, x') 



m e ao 



S(x — x') (hardcore) 
1 

(Coulomb) 



(19) 



\x — x'\ 



As shown in Figure [5] and anticipated before, the sharp- 
ening of the Wannier function at the origin causes greater 
interaction matrix elements, when pressure is applied, ir- 
respective of the chosen electron-electron interaction. 

As previously stated, the on-site electron-electron in- 
teraction V from above equals the Hubbard U, since 
higher energy orbitals, and thus screening effects, are ab- 
sent by construction. 
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FIG. 4: (Color online) Nearest neighbor transfer integral t = 
{ipR = o\Ho\tpH^ a ) for a model with Dirac-comb potential, as a 
function of lattice constant a/an. 



FIG. 5: (Color online) On-site electron-electron interaction 
V, for a hardcore interaction (main graph), and a Coulombic 
one (inset), as a function of lattice constant a/an. 
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FIG. 6: (Color online) Wannier function Vo(^) and corre- 
sponding on-site hardcore interaction (see Eq. (|18[0 as a func- 
tion of the gauge parameter c (see text for definition). A 
maximally localized Wannier function is obtained for c = 0.0 
(solid lines) and c = 1.0 (dotted). For c = 0.5 (dashed), 
the Wannier function is symmetric with respect to a/2, and 
has equal weight at the site R = and R = a. The lattice 
constant is fixed at a/ao = 5.0. 



6. Maximally localized Wannier functions versus maximal 
Hubbard interaction 

In Sec. I we mentioned that another technique to 
choose the Wannier function gauge is given by the request 
to maximize the static local, partially screened Coulomb 
interactio n 19 ! 25 - the Hubbard U. 

While not actually performing this approach for our 
model, we give evidence that in this simple case, both 
techniques are equivalent. As said before, the Hubbard U 
equals the bare Coulomb interaction V since we consider 
only a single band, so, contrary to the general casei 9 -, the 
argument does not involve any screening related effects. 

The quest is thus to find a Wannier gauge, meaning 
an additional factor exp(z</>fc) in Eq. (|10p . that yields the 
greatest interaction element as given by Eq. (|18p . The 
choice of gauge can be absorbed into the distribution A ni 
and we define 



M<t>\ = 



v/a dk 



C[4>] J-ir/a 2^ y/l + 2^ =1 s m co8(kma) 

(20) 



Having seen the correspondence between the value of 
Aq and the interaction strength, one might endeavor to 
solve the functional derivative 5Ao[0]/i50 = for </>. Yet 
the dependence of the normalization constant, C, given 
in Eq. (fT^ . makes this a rather tedious task analyti- 
cally. Instead, we shall argue for a specific example for 
the choice of Wannier gauge, and make some general 
comments. Consider a gauge field that is linear in mo- 
mentum, 4>k = —cka [see also Refill]. Figure [5] displays 



the Wannier function (upper panel) and the resulting on- 
site interaction (bottom panel) for the hardcore case (see 
Eq. (fl9|l ). for different gauge parameters c, and for a fixed 
lattice spacing a/ao = 5.0. With c =/= the inversion sym- 
metry of the Wannier function with respect to the site 
R = is lost. This was the requirement for maximally 
localized Wannier functions in one dimensio n 14 ' 24 , and as 
seen in Figure for c > the tail of the Wannier func- 
tion augments, and the Coulomb interaction decreases. 
While c = 1 is a mere lattice translation, c = | amounts 
to a case, where the Wannier function of site R = has 
equal positive weight at x = and x = a, corresponding 
to A® = A\. As discussed above, an increase in Aq is 
caused by the mixing in of negative components to tpo 
from neighboring sites, leading to a smaller normaliza- 
tion factor C in Eq. (|14[) . Owing to the symmetry, in 
the current example, c = ^, negative contributions A m 
will come only from the sites 2m + 1 for m < and 
2m for m > 0. Yet the overall gain in renormalization 
is distributed over the two equivalent positions m = 0, 
and 1 for which A m > 0. As a result of this shifting 
of weight to the site m = 1, Aq and thus the on-site 
Coulomb interaction V are much lower than in the case 
with inversion symmetry around that site. As seen in 
the bottom panel of Figure O it is indeed the maximally 
localized Wannier functions (c = 0.0,1.0) that yield the 
greatest possible Coulomb interaction matrix element for 
our simple model. 

For gauge fields <f>k = —ck a a depending on the momen- 
tum to the power a > 1, the argument is geometrically 
less obvious, but still true as verified numerically. Indeed, 
only the integrand in the coefficient A n= o of Eq. (f^Uj) is 
always positive for <pk — as a function of k. Therefore 
any modulation in cos(0fc), with <fik ^ 0, will decrease 
the corresponding integral to a greater extent than for 
n 0, in which case the integrand changes sign with k 
already for <pf~ = 0. As a consequence, the decrease in the 
n = contribution to the overall normalization C will be 
greater on a relative scale than for n / 0, and thus Aq 
decreases with any non-constant <f>k ^ 0. 

This can be taken as a further indication that also for 
realistic systems, the maximally localized Wannier func- 
tions and the maximal Hubbard U approach are generally 
giving the same results. 



B. fee hydrogen under pressure 

Towards a more realistic application of the gained in- 
sight, we apply the fully ab initio approach^ 9 - of the 
cRPA within maximally localized Wannier functions to 
the "simplest" realistic system, namely solid hydrogen. 

While at low pressure, solid hydrogen forms a molecu- 
lar crystal. It was conjectured, already in the 1930s, that 
at high pressure hydrogen should become an isotropic 
metals. Here, however, we shall not be interested with 
the precise phase diagram of solid hydrogen. For the 
sake of simplicity, we assume throughout the discussion 
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FIG. 7: (Color online) LDA band-structure of fee hydrogen for 
the two lattice spacings a = 5.0 a.u. (solid) and a = 10.0 a.u. 
(dashed) . The Fermi level corresponds to the origin of energy. 



a face-centered cubic crystal structure with variable lat- 
tice constant. Moreover, we are well aware that for the 
current case the problem of self-interaction22 within the 
LDA formalism is a particularly severe issue. However, 
here, we are not concerned with the accuracy of the LDA 
band-structure but with compression induced trends in 
a 3d multi-orbital setup. 

For the one-particle band-structure, we employ in this 
work the LDA in the full-potential FP-LMTO^ realiza- 
tion. In the LMTO basis, we include orbitals up to the 4f, 
using local orbitals^ for multiple orbitals per 1-channel, 
and use a Brillouin zone discretization with up to 10 3 
points. As described in Ref^, maximally localized Wan- 
nier functions are then constructed for the hydrogen Is 
band, which is entirely isolated from all other Kohn-Sham 
excitations. In other words, the effective model consists 
of a single half-filled orbital. We stress that while the 
Wannier functions of different sites are orthogonal by 
constructions, the LMTO basis functions - in analogy to 
the tight-binding paramctrization in the preceding sec- 
tion - are not. As a consequence, the same prototypical 
response to pressure as discussed above can be expected 
in the current case. 

Figure [7] shows the LDA band-structure for the ex- 
tremal lattice constants that we consider. As expected, 
under growing compression, the dispersions increase and 
unoccupied bands are shifted upwards 3 ^. 

In line with this is the behavior of the hopping t — the 
nearest neighbor transfer matrix element in basis of the 
maximally localized Wannier function for the subspace of 
the lowest Kohn-Sham excitation. As shown in Figure [51 
it augments with decreasing lattice constant a, account- 
ing for the greater derealization. In Figure [8] is shown 
the maximally localized Wannier function (it is real, cf. 
RefjH) of the hydrogen atom at the origin, as a function 
of the (scaled) distance £ towards the nearest neighbor at 



e.g. r = (£ = 1,0, £ = l)a/2. When the lattice constant 
shrinks, clearly witnessed is both, a growing weight at 
the neighboring position (£ = 1/2), as well as an increase 
in the Wannier function at the origin (£ = 0). This is in 
complete analogy to what was shown in Figure [2] for the 
model considered above^ 3 .. 

It is thus expected that the on-site matrix element of 
the bare Coulomb interaction Vr,r grows under com- 
pression as before. And, indeed this is the case, as can 
be inferred from Figure [5] in the second panel from the 
bottom. Yet, does that entail for the Wannier spread f2 
the same behavior as witnessed in the one dimensional 
model ? In three dimensions, the gain in spread by ef- 
fects of hybridizations with neighboring sites might turn 
out less prominent, since the region occupied by nearest 
neighbor atoms is relatively small. Hence, the angular 
integral in O, Eq. ([I]), even for the radius of the distance 
to the 12 neighboring atoms (fee) will run very much over 
a sphere on which the Wannier function is mostly zero. 
And, indeed, as displayed in the second panel from the 
top in Figure [HI the Wannier spread actually decreases 
under compression 3 ^. 

A direct interpretation of the spread in pressure stud- 
ies has thus to be taken with caution. First, depending 
on the crystal-structure (dimension, number and position 
of neighboring atoms) the spread does not necessarily re- 
flect the increased derealization of charge carrier but 
can be dominated by the more isotropic concentration of 
weight at the origin (sec Figure [5]) . Second, as seen in 
Figure [9] (and already in the seminal work, Ref Table 
I), the momentum space convergence of f2 is poor. Also, 
a change in lattice constant upon compression, changes 
the accuracy of the spread function when the number of 
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I with r=(£,0£)a 

FIG. 8: (Color online) Maximally localized Wannier function 
of the origin, as a function of the scaled distance f towards 
the nearest neighbor atom at r/a = (£ = i,0, £ = |). The 
arrow in the upper left indicates lf-^/n, which is the value of 
V>(0) for the exact atomic Is hydrogen orbital (see text for 
comments). The inset shows a zoom on the position of the 
nearest neighbor atom. 10 3 k points were used. 
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FIG. 9: (Color online) Maximally localized Wannier functions 
and cRPA for the Is orbital of fee hydrogen. Shown are (from 
top to bottom) the nearest neighbor transfer integral ("hop- 
ping") t, the spread fi of the Wannier function (see Eq. @), 
the bare on-site Coulomb interaction Vr,r, and the cRPA 
screened on-site "Hubbard" interaction U as a function of the 
lattice constant a, and for different discretization of the Bril- 
louin zone (4 3 , 8 3 , 10 3 , and for the spread : 20 3 ) as indicated. 
The symbols in the screened interaction U indicate the val- 
ues obtained when substituting in the cRPA polarization the 
LDA eigenvalues with those of the calculation for a/eto = 5.0 
(using 10 3 k-points). 



k-points is kept constant. While the spread is of course 
the entity that is minimized in order to obtain the max- 
imally localized Wannier functions for a given pressure, 
the spread itself is not a reliable measure for trends in 
the Wannier functions upon compression. Instead, one 
should either plot the functions, or resort to the bare 
Coulomb interaction matrix elements, which - owing to 
the different powers of the position operator - converges 
well for a moderate k-point sampling. 

In the current case, while still constructing an effec- 
tive oTie-band model, the initial system contains higher 
energy orbitals. As a result, and contrary to the simple 
tight-binding model, there is a non-zero partial polar- 
ization, P r of Eq. that screens the bare Coulomb 
interaction V. The on-site part of the screened interac- 
tion Wtq o - the Hubbard U - within cRPA is displayed 
in the lowest panel of Figure [5] 

Interestingly, this quantity, in contrast to the bare in- 
teraction, shows a non-monotonic behavior, that reflects 
the struggle of two opposing effects in the polarization 



P r . As can be inferred from Eq. ([2]), and the equation 
for P s , Eq. ([3]), the changes in the polarization P r origi- 
nate from modifications of transitions from the Wannier 
"Is" into orbitals at higher energy. These can be altered 
by two ingredients (see again Eq. ([3])) : the transition 
matrix elements, and the (Kohn-Sham) transition ener- 
gies ekn, i-e. the band-structure. The effect of pressure 
will be different for these two mechanisms. Indeed, the 
increasing compression pushes the bands further apart, 
as seen in Figure El diminishing the polarization. The in- 
creased overlaps/hybridizations of orbitals, on the other 
hand, tend to make the polarization grow. In order to 
separate the influence of the two contributions, we com- 
puted the partially screened Coulomb interaction, U , as a 
function of lattice spacing, albeit while keeping the band- 
structure fixed at the values obtained for the highest pres- 
sure, a/ao = 5.0. Thus the corresponding screened inter- 
action will contain only the effect of changes in the Bloch 
functions. As indicated by the symbols in the figure of 
the Hubbard U, the latter decreases under compression, 
proving the above conjectured opposition in trend to the 
influence of the band-structure. 

This effect is surely very material specific. In systems 
with more electrons, where pressure will e.g. cause the 
enhancement of bonding/anti-bonding splittings, it can 
be expected that the changes in the band-structure often 
prevail such as to reduce the polarization under compres- 
sion. 



IV. CONCLUSIONS AND PERSPECTIVES 

In conclusion, we have studied the influence of exter- 
nal pressure onto the construction of effective low energy 
many-body systems. Using maximally localized Wannier 
functions for a one-band tight binding model, we ratio- 
nalized the counter-intuitive, yet prototypical behavior 
of the bare Coulomb interaction, namely that its matrix 
elements augment upon compression, as a consequence 
of the derealization of the Wannier functions. This 
we understood to be caused by increased admixtures of 
non-orthogonal nearest neighbor tight-binding functions 
when the lattice spacing shrinks. As a more realistic sys- 
tem, we investigated fee hydrogen under pressure, and 
constructed an effective model for the half-filled Is or- 
bital using maximally localized Wannier functions. For 
the transfer integrals and the bare Coulomb interaction, 
we witnessed the same tendencies as in the model case. 
Yet, the Hubbard U, calculated as the on-site screened 
interaction within the constraint RPA technique, exhib- 
ited a non-monotonous trend under compression. This 
we traced back to a struggle between two opposing ef- 
fects in the strength of screening. All these highlight the 
intricacy of mechanisms that influence effective models, 
and emphasizes the need for reliable ab initio techniques 
for their construction. 
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The fact that the energy difference between the Is and 
2s, 2p bands at a — lO.Oa.u. is already smaller than the 
atomic ionisation energy is owing to the aforementioned 
spurious self-interaction in the LDA. 

We note that the issue of self-interaction inherent in the 
LDA treatment impedes the Wannier function to converge 
towards the well-known Is atomic eigenfunction ip(r) — 



\ / y/ ndQexp(— r/ao). This is why the values of il>(0) for 
large separation are below teh indicated 
We note that the true value of the spread at infinite sepa- 
ration should be 3oq - the atomic value. Again this is due 
to the inadequateness of the LDA. 



